A Combined Quantum Mechanical and Molecular Mechanical Potential for Molecular Dynamics Simulations

A combined quantum mechanical (QM) and molecular mechanical (MM) potential has been developed for the study of reactions in condensed phases. For the quantum mechanical calculations semiempirical methods of the MNDO and AM1 type are used, while the molecular mechanics part is treated with the CHARMM force field. Specific prescriptions are given for the interactions between the QM and MM portions of the system; cases in which the QM and MM methodology is applied to parts of the same molecule or to different molecules are considered. The details of the method and a range of test calculations, including comparisons with ab initio and experimental results, are given. It is found that in many cases satisfactory results are obtained. However, there are limitations to this type of approach, some of which arise from the AM1 or MNDO methods themselves and others from the present QM/MM implementation. This suggests that it is important to test the applicability of the method to each particular case prior to its use. Possible areas of improvement in the methodology are discussed.